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IMPINGEMENT OF WATER DROPLETS ON WEDGES AND DIAMOND 
w AIRFOILS AT SUPERSONIC SPEEDS 

03 

^ By John S. Serafini 


SUMMARY 

An analytical solution has heen obtained for the equations of motion 
of water droplets in^jinging on a wedge in a two-dimensional supersonic 
flow field with a shock wave attached to the wedge. The closed- form 
solution yields analytical expressions for the equation of the droplet 
trajectoiy, the local rate of in^jingement and the in5>ingement velocity at 
any point on the wedge surface, and the total rate of impingement. The 
a na lytical expressions are utilized in the determination of the impinge- 
ment of water droplets on the forward surfaces of diamond airfoils in 
supersonic flow fields with attached shock waves. 

For a wedge, the lesxilts provide Information on the effect of the 
semi-apex angle, free- stream Mach number, pressure altitude, and droplet 
size. For the diamond airfoil, additional calculations provide informa- 
tion on the effect of airfoil thickness ratio, chord length, and angle 
of attack. 

The results for the diamond airfoils are also correlated in terms of 
the total collection efficiency as a function of a modified inertia 
parameter. The results are presented for the following range of vari- 
ables: droplet diameter from 2 to 100 microns, free- stream Mach n-umbers 

from 1.1 to 2.0, pressure altitudes from sea level to 30,000 feet, free- 
stream static temperatures from 420° to 460° R, semi-apex angles for the 
wedge from 1.14° to 7.97° and corresponding diamond airfoil thickness-to- 
chord ratios from 0.02 to 0.14, and chord lengths from 1 to 20 feet. 


INTRODUCTION 

The problem of ice prevention on aircraft flying at subsonic speeds 
up to flight critical Mach numbers has been a subject of considerable 
study and research by the NACA. The recent advent of aircraft flying at 
transonic and supersonic speeds has required an extension of these icing 
studies to higher speeds. That an icing problem exists in the transonic 
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and supersonic speed range is verified in reference 1, wdich shows hy an 
analytical investigation with experimental confirmation that diamond 
or symmetrical double-wedge airfoils are subject to possible icing at 
flight Mch numbers as high as 1.4. A similar result is expected for 
other airfoil shapes that are being considered for use at treinsonic and 
supersonic flight speeds. 

In conducting research on the problem of ice prevention on aircraft 
and missiles, regardless of the magnitude of the flight speed, it is 
essential that the impingement of cloud droplets on airfoils and other 
aerodynamic bodies be determined either through theoretical calculations 
or through experimentation. The in^jingement variables which must be 
determined are the total water catch, the extent of impingement, and the 
rate of in 5 )ingement per tmit area of body sxarface. These variables can 
be determined analytically from calculations of the cloud-droplet trajec- 
tories obtained for the various aerodynamic bodies. Investigators have 
reported the results of cloud-droplet trajectories about right- circular 
cylinders (refs. 2 to 5) and about airfoils (refs. 6 to 9) immersed in 
an incompressible fluid. An evaltiation of the effect of con^pressibility 
on the droplet trajectories about cylinders and airfoils up to the criti- 
cal flight Mach number is presented in reference 10. 

At present there exists little information on the impingement of 
droplets on aerodynamic bodies in a supersonic air stream. Concentration 
of pant effort on the problems of inpingement on airfoils at subsonic 
flight speeds and the present lack of convenient and rapid means for 
obtaining the rotational flow fields about airfoils at supersonic speeds 
are possible explanations for the scarcity of trajectory calculations for 
the supersonic region. An initial contribution to the solution of the 
over-all problem of impingement of water droplets on aerodynamic bodies 
at supersonic speeds is given in reference 11, which presents an analysis 
of the water- interception characteristics of a wedge in a supersonic flow 
field. 

The present report extends the analysis of reference 11 and further 
presents an extensive study of the inpingement of water droplets on two- 
dimensional wedges and diamond airfoils for sipersonic flight speeds that 
result in attached shock waves and constant velocity fields behind the 
shock waves. For the wedge angles and diamond airfoil thickness ratios 
to be considered herein, the shock-wave attachment Mach number varies 
from a value slightly greater than 1 to aboTxt 1.4. The method enployed 
is based on an analytical solution of the equations of motion by means of 
a closed-form integration. The closed-form solution yields analytical 
expressions for the equation of the trajectories, the local inpingement 
efficiencies, the velocity at any point on the trajectories, and the total 
rate of inpingement. This solution has been made possible by using an 
enpirical relation for the drag coefficients for spheres that gives a 
good approximation to the exp^imental drag coefficients. 
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The results of calculations for the rate^ extent, distribution 
of the impingement of water droplets on wedges and diamond airfoils are 
‘ presented herein. The ranges of variables incliaded for the wedge are 

Mach number from 1.1 to 2.0; pressure altitude of sea level, 15,000 feet, 
and 30,000 feet; droplet diameter from 2 to 100 microns; free- stream 
static ten 5 )erature of 420°, 440°, and 460° R; and the tangent of the semi- 
apex angle from 0.02 to 0.14. The ranges of variables for the i nmnnfl or 
double-wedge airfoil are the same as those for the wedge, and the addi- 
co tlonal variables for the diamond airfoil range from 1 to 20 feet for the 

m chord length and from 2 to 14 percent for the thickness . 

The work presented in this report was performed at the MCA Lewis 
laboratory. 


AMLYSIS 

Statement of Problem 

The solution of the problem of Impingement of water droplets on a 
two-dimensional wedge at supersonic speeds with an attached shock 
wave is not as difficult as that for the inipingement on various airfoils 
at low subsonic speeds . For the wedge at supersonic speeds with an 
attached shock wave, the air velocity everywhere ahead of the shock wave 
is constant and equal to the free-stream air velocity Vq_ (fig. l). The 

air velocity behind the shock wave is also everywhere constant and 

parallel to the wedge surface. All the droplets have the same initial 
velocity (that of the free-stream air velocity), and their trajectories 
are exactly coincident with the air streamlines upstream of the shock 
wave. All water droplets of a given size are subjected to identical air- 
velocity fields, which in turn produce identical force systems downstream 
of the shock wave. Irrespective of the point along the shock wave where 
the droplets cross the wave. It follows, therefore, that, for droplets 
of a given size, all the trajectories in a given problem are identical 
with respect to the point where the droplet crosses the shock wave. 

By adopting a frame of reference which moves at the constant velocity 
of the air Vg downstream of the shock wave, the problem of the droplet 

motion is reduced to the still- air problem, defined as the determination 
of the motion of a droplet that, having an initial velocity, is projected 
into qiilescent air. Hence, relative to the moving frame of reference, the 
initial velocity of the droplet upon crossing the shock wave is equal to 
the vectorial difference of the free-stream air velocity the air 

velocity Vg downstream of the shock wave. Adoption of the frame of 

reference moving with a constant velocity reduces the problem from the 
solution of two simultaneous nonlinear second-order differential equations 
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in the fixed coordinate system to the solution of a single nonlinear 
second-order differential equation in the moving coordinate system. At 
any given instant, the droplet displacement relative to the point of 
intersection with the shock wave in the fixed frame of reference is 
obtained hy adding vectorially the droplet displacement within the moving 
frame of reference to the displacement of the moving frame of reference 
for the same increment of time. 

This general method of attack was used in reference 11, where the 
one second-order differential equation representing the droplet motion 
relative to the air velocity behind the shock wave was integrated graphi- 
cally. However, it is possible to obtain a coinpletely analytical solution 
by means of a closed-form integration without resorting to the use of 
numerical integrations or analog conputing equipment, if an empirical 
relation is assumed for the drag coefficient as a function of the Reynolds 
nimiber of the droplet relative to the air. It will be shown that this 
closed-form integration of the still-air problem when applied to the wedge 
in supersonic flow with attached shock wave yields the equations for the 
trajectories of the water droplets and the droplet velocities at any point 
on the trajectories and makes available relations for the rates of total 
water inpingement and the local rates of water impingement along the wedge 
surface. Furthermore, it is shown that these equations can also be 
readily applied to the determination of the droplet Impingement on a 
diamond airfoil in supersonic flow with attached shock waves. 

Most of the -usiial assumptions made in the previous investigations on 
impingement at subsonic speeds are also required for this investigation. 
These assunptions are: (l) The water droplets are always spherical and 

do not change in size, (2) the force of gravity on the droplet may be 
neglected in conparison to the drag forces, (3) the drag of the air on 
the droplet is that of a viscous inconpressible fluiid. Here it is addi- 
tionally assumed that (a) the two-dimensional supersonic flow field about 
the wedge is frictionless except within the infinitesimally thin attached 
shock wave, (b) no condensation shock occurs and no change in phase occurs 
as the water droplets traverse the oblique shock wave, and (c) the 
unbaleince of the forces on the water droplet from the instant it enters 
the shock wave until it emerges from the shock wave can be neglected in 
the calculation of the trajectories. 


Equation of Droplet Mation in Moving Reference Frame 
The velocity of the droplet in the moving frame of reference is 



U = 


( 1 ) 
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where is the droplet velocity with respect to the fixed frame of 

reference, and Vg is the air velocity downstream of the attached shock 

wave also with respect to the fixed frame of reference (fig. l). In the 
frame of reference moving with the velocity V2 (air velocity downstream 

of shock wave), the statement of Newton's law of motion for the water 
droplet becomes 

D = - Cd (I) P2«aV = (I) «a^p^ 

from which 

A complete list of symbols is given in appendix j 

Equation (3) is the differential form of the equation of motion of 
a droplet which is projected with an initial velocity into a region of 
quiescent air (the so-called still-air problem). The shock wave is con- 
sidered to be a surface of discontinuity from which the droplets emerge 
with a velocity V2^. In this case the initial velocity of the droplet is 




which is the magnitude of the vector difference of the air-velocity vectors 
upstream and downstream of the attached shock wave. As can be shown from 
a consideration of the continuity equation and the equation for conse^a- 
tlon of momentum across the oblique shock wave, the velocity vector 

is normal to the shock wave. At a^ subsequent instant of time the 
relative droplet velocity vector U retains the same angiilar orientation 
to the shock wave and changes only in magnitude. 

In reference 11 the solution of equation (3) has been obtained by 
numerical integration. The result obtained in this manner makes it neces- 
sary to use a graphicel procedure in determining the trajectories and the 
local rates of impingement. However, an analytical solution of equa- 
tion (2), which eliminates the graphical procedure, can be obtained if 
the experimental values of the drag coefficient Cj) are expressed in a 

function involving the Reynolds niimber Re^. The relation is 
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Cd =||: (1 (5) 

where € and m are the empirical constants. This empirical relation 
is a valid approximation in the range of Reynolds numhers to which cloud 
droplets are subjected in trajectory calculations. Substitution of the 
expression for C-q (eq. (s)) in equation (3) results in the expression 



where the local relative Reynolds number Re^ = 2p2Ua/ii2- The displace- 
ment of the water droplet in the moving frame of reference is x and is 
measured from the air streamline which intersects the shock wave at the 
point where the water' droplet entered the air-flow field downstream of 
the shock wave. The closed-form integration of the differen'tial equa- 
tion (6) is presented in appendix B. The use of 2/3 for the exponent 
m and 0.158 for the value of the empirical constant e in eqiaation (5) 
yields an empirical curve for the drag coefficient as a function of the 
local Reynolds numbers that approximates very well the variation of the 
experimental values of the drag coefficient in the range of Reynolds 
numbers from 0.5 to 500. The value of 2/3 for the exponent m also 
facilitates the closed-form Integration of the differential equation of 
motion. In figure 2 a graph of the empirical relation is presented, 
along with the drag coefficient data of references 4 and 12. 

The results of the integration are given by the following eq-uatlons: 


- tan"^ 


(V) 

where 

T dimensionless time variable, (bp^/p^a^)! 
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initial value of the local relative Reynolds number Re^.^ 

2ap2Ui/ \i2 

a droplet radius, ft 

water density, 1.9398 slugs/cu ft 

pg density of air behind shock wave, slug/cu ft 

e empirical constant of relation of drag coefficient as function 

of Reynolds number, 0.158 


and 


where 








cp = tan"^ 


0 ^ cp ^ rt /2 


The intermediate steps of integration are given in appendix B. 


( 8 ) 


( 9 ) 


Equation (7) and (8) give, respectively, the displacement and the 
velocity of the droplet at any instant in the moving frame of reference. 
The displacement of the droplet with respect to the point where it crossed 
the shock wave can be obtained by a vectorial addition of the displacement 
X and the displacement of the moving reference frame in the corresponding 
time interval. The droplet velocity relative to the fixed frame of 

reference must ^so be obtained by the vectorial addition of U 
(eq. (8)) and Vg (the air velocity downstream of the shock wave). Equa- 
tion (9) gives the maximum value of x that is obtained as the time of 
travel in the air-flow field downstream of the shock wave approaches 
infinity. The significance of this qtiantity will be discussed in sub- 
sequent sections. 
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Relations Required for Application of Closed-Form Solution 
to Obtain Droplet Motion and Inpingement in 
Fixed Reference Frame 

Impingement on vedges . - For a problem of given aerodynami'^ condi- 
tions, tbe trajectories of al T the water droplets for a given size are 
identical when the points where the droplet trajectories intersect the 
shock wave are superimposed. This unique characteristic of the water- 
droplet trajectories about a wedge in supersonic flow with an. attached 
oblique shock wave is the result of two constant velocity fields, one 
upstream and one downstream of the shock wave. Therefore, only one set 
of equations for a single trajectory will be necessary to calculate the 
inpingement parameters for a specified problem, including a given droplet 
size. The values of the initial relative velocity the initial 

Reynolds number Re^ and the density ratio needed for 

substitution in the closed-form solution of the equations of droplet 
motion. These values can be obtained from information available in 
reference 13 and from the use of simple algebraic and trigonometric rela- 
tions for given v^ues of the free-stream Mach number the free- stream 

static tenperature t^, the angle of siarface inclination to the free- 
stream direction o, the free-stream static pressin:e pj_, and the droplet 

diameter d. These relations result in the following expressions for the 
initial relative velocity and initial relative Reynolds number; 

Ui = I U^l = -Jl +Q2 _ 2Q cos a = (lO) 

Rer,i = 2 

where 
Q Vg/Vi 

R 53.3 ft-lb/(lb)(OF) 

IT 1.4 

g 32.2 ft/sec2 
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A convenient form of the solution for the impingement on a wedge or 
the front half of a diamond airfoil is obtained if S is defined as the 
distaince to the point of impingement measured from "Che leading edge for 
a water droplet that enters the flow field behind the shock wave at a 
distance ^ above the leading edge (fig- 1). The unique relation between 
S and ^ in a given problem for droplets of the same size is quite 
readily determined by considering the displacement of the water droplets 
as the vectorial sum of the displacement of the water droplet relative 
to the moving frame of reference and the displacement of the moving refer- 
ence frame relative to a fixed frame of reference (referred to wedge). 
Since the moving reference frame has a velocity equal to the air velocity 
V 2 , which is constant in magnitude and parallel to the wedge surface, 

only the droplet travel in the moving reference frame contains the com- 
ponent of droplet travel which represents the approach of the water 
droplet to the wedge surface. For a water droplet starting from point A 
and impinging on the wedge s\arface at point D (fig. l), the displacement 
of the moving reference frame (= V 2 t, where t is zero at point A) is 

given by the displacement vector AB equal to CT, and the droplet motion 
in the moving frame of reference is given by the displacement vector M 
equal to AC. Therefore, relative to the starting point at A (fig. l), 
the displacement of the water droplet to the point of impingement at D 
is obviously eqtial to (AB + BD) or to (AC + .CD). From figure 1 the 
displacement of the water droplet at the point of impingement D, meas- 
\ired from the leading edge at E, is given by adding the vector EA to 
the displacement vector from the starting point A, and the displacement 
of the droplet at D referred to the leading edge is 


s = |ed| = |ea + AC + cd| = [eg + cd| = |ec| + |cd| = i + i' (12) 

where ^ is defined as the magnitude of the displacement vector K and 
(= V 2 t) is defined as the magnitude of the displacement vector 

(the displacement of the moving frame of reference). As previously noted, 
S is the distance from the leading edge to the point of inpingement of a 
water droplet that initially (upstream of the shock wave) was on the 
streamline at a distance ^ from the leading edge, measured normal to 
the streamlines. 

The values of ^ and i are obtained in terms of x, the distance 
of travel in the moving reference frame, from simple trigonometric identi- 
ties involving the various angles shown in figure 1: 


^ = X sin 9 tan (v+cr) (13) 

where x is given by equation (7), and v is the angle between the rela- 
tive velocity vector U and the free-stream air velocity vector and 

is given by 
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(13a) 


I = X sec (V-KJ) (l4) 

Substitution of equation (14) into equation (l2) for the s\irface distance 
to the point of inrpingement yields the following: 

S = X sec (i> - ht) + V 2 (15) 

where t is the dimensionless time variable used in the determination 
of X and is defined in the previous section. 

Since x is a function of T in equations (13) and (ib)^ the expres- 
sions for ^ ani^ S, respectively, are functions of T. However, since 
T cannot be eliminated from equations (13) and (15), S cannot be obtained 
explicitly as a function of C- Nevertheless, the curves of the 
initial displacement of the -trater-droplet trajectory from the leading edge 
normal to the free-stream direction, against S, the distance to the point 
of ingpingement of the stated water-droplet trajectory, can be obtained by 
substitution of the same set of values for T in the expressions for ^ 
n.nfl S. 


An analytical expression for the local impingement efficiency, which 
is designated an P, can be obtained from the above expressions for ^ 
and S. The local in^ingement efficiency P is defined by the expression 


P = 




dT 

dT 


(16) 


where ^ is the difference in the Initial displacements of two water 
droplets having very nearly eq^uaJ. initial displacements, and ^ is the 
small increment of wedge siirface between the points of impingement of the 
two water droplets- Differentiating ^ and S with respect to T and 
performing the division indicated by eq\aation (l6) yield the following 
expression for P: 


°2 ^ ^ 6eV^ ^ 
ni ni a fv 




-2/3 


+ 1 1 e' 


1 


zjz 


(17) 
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where 


iijL = sin 9 tan (v+a) (l8a) 


H 2 = sec (h+o) (l8b) 

p a^ 

“3 ^ ^2 3 ^“ = VgCt/r) (l8c) 

Since p and S both are functions of T, the local impingement effi- 
ciency P at any point on the wedge surface is determined by using the 
same value of T in equations (17) and (l5). The value of p that 
exists as the point of impingement of the water droplet on the wedge sur- 
face approaches the leading edge as a limit (S— »-0) is defined and 
given by the following: 


Po = 



sin a 


(19) 


The magnittide and direction of the droplet velocity at the point of 
in^jingement (relative to the fixed frame of reference) can also 

be easily obtained at any point on the wedge surface as 


^d,im - ^I'l 


+ 0)2 + 2a CD cos (v+a) 


( 20 ) 


•^m = ^ - crijn = ^ 


sin"^ 



sin (u+a) 


( 21 ) 


where K is the eingle between the free- stream velocity vector nnd 

the droplet velocity vector Vd. and ex' is the angle between the droplet 

velocity vector and the air velocity vector Vg. In equations (20) 

and (21) for a given trajectory, ^2, ox, V, and a are constants. 

Therefore, for a given trajectory, and are fxinctlons only of 

U (eq. (8)) -wbicb is in turn a function of t, the dimensionless time 
variable. 



12 


KACA TN 2971 


Impingement on diamond alrfoila . - The In^lngement on a diamond air- 
foil may be obtained from the solution to tbe problem of impingement on 
a wedge as presented heretofore. In this report a diamond airfoil is 
considered to be a symmetrical double-wedge airfoil, the maximum thickness 
of which occurs at 50 percent of chord (fig. 3). At zero angle of attack, 
the impingement on a diamond airfoil will be limited to the region from 
the leading edge to the shoulder at 50 percent of chord. The solution 
for impingement on a wedge surface having a given semi- apex angle a (the 
angle of inclination of either wedge stnrface to the free- stream direction) 
can also be used as the solution for a diamond airfoil where the thickness 
ratio of the diamond airfoil is equal to tan 0 , the tangent of the semi- 
apex angle, and where the droplet size and other parameters of the problem 
are the same as for the wedge. Therefore, the values of the local impinge- 
ment efficiencies 3 and at any given point on the surface will be 

identical for both the wedge and diamond airfoil under the aforementioned 
similarity of conditions. 


The solution for the Impingement on the diamond airfoil at angle of 
attack can also be obtained from the solution for Impingement on a wedge 
as for the case of the diamond airfoil at zero angle of attack. l-Jhen 
the symmetrical diamond airfoil is at angle of attack a, the angle of 
inclination of its forward ig)per sinrface to the free- stream direction is 
equal to (a-a) and that of the forward lower 'surface of the diamond 
airfoil is equal to (a+a). Therefore, the solution to the impingement 
on the upper and lower surfaces of the diamind airfoil is obtained from 
the solutions for impingement on wedges having the redefined semi- apex 
angles of ( 0 -a) and ( 0 +a), respectively, where the droplet size and 
other parameters of the problem are kept the same. For the diamond air- 
foils at angles of attack having tangents equal to or greater than the 
thickness ratio, the water droplets will not Impinge on the upper surface. 
At angles of attack having tangents greater than the thickness ratio, 
some water droplets may impinge on the lower surface beyond 50 percent of 
chord. These conditions are presented in figure 3, which illustrates 
schematically the three conditions. For a< tan"^ (t/c), the Impinge- 
ment occurs on surfaces AC and AB; for a = tan“l (t/c), impingement 
occurs only on sirrface AC; for a> tan“^ {t/c), in^1 ngement occurs on 
lower surface AC and may occur on lower surface CD. However, the 
condition where a> tan”l (t/c) is not considered herein, since the 
solution presented in this report is not valid for the determination of 
trajectories of droplets impinging on the surface beyond the shoulder or 
50 percent of chord of the diamond airfoil (sturfaces BD or CD) where 
a portion of the trajectories is within the expansion zone emanating from 
the shoiolder. 
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RESULTS AMD DISCUSSION 

From ttie equations presented in the previous section and in appen- 
dix the iHg)ingement of water droplets on the wedge in a supersonic 
flow field with an attached shock wave can he calculated over a range of 
free-stream conditions, wedge angles, and droplet sizes. As has been 
previously indicated, the impingement characteristics of the diamond 
airfoil (at zero angle of attack and also for small values of angle of 
attack) can be readily determined from the in5>ingement on wedges for 
similar conditions. The results for the wedge and diamond airfoil are 
presented and discussed separately. A comparison of the total collection 
efficiency and the water ingoingement rate for a diamond airfoil at zero 
angle of attack with those for a HACA 0006-64 airfoil at zero angle of 
attack is presented in appendix C. 


Wedge 

Local impingement efficiency . - The rate of water impingement on a 
local area of wedge or airfoil surface is proportional to a dimensionless 
term P, the local impingement efficiency. The local rate of water 
impingement in pounds per hour per square foot is 


Wp = 0.329 T]_wp 

where p is the fractional part of the maximniiii water that could impinge 
on a local area of the wedge or airfoil, if all the droplet trajectories 
were parallel to each other and the surface of impingement were projected 
into a plane perpendicular to the trajectories. The local impingement, 
when given as a function of the surface distance of the wedge, allows the 
determination of the local rate of impingement of water droplets at any 
point on the surface, the total impingement of water droplets on the 
entire surface or any given portion of the surface, and the extent of 
impingement on the siirface. The local impingement efficiency P is 
related to a point at a given distance S on the wedge surface in eq\ia- 
tions (l5), (17), and (l9) by the dimensionless time variable which 
is common to 1 three expressions. The variation of the local impinge- 
ment efficiency P with S, the s\irface distance along the wedge, is 
presented in figure 4 for an extensive range of free-stream conditions, 
semi-apex angles, and droplet sizes. The values of p at the leading 
edge (S— >0) is the sine of the semi-apex angle (sin a); and as S 
increases, the value of P decreases rapidly and approaches the value of 
zero asymptotically as S approaches Infinity. However, it is to be 
noted that negligibly small values of p (p « 1 percent of Pq for 
the wedge) are attained at large but finite values of S. 
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The curves of P as a function of S presented in figinre 4 are 
those of an idealized situation. The assumed two-dimensional supersonic 
flow field about the wedge does not account for a stagnation point that 
must exist at the leading edge of the wedge, regardless of the sharpness 
of the leading edge. In addition, the leading edge of wedges and diamond 
airfoils might he considered to he somewhat rounded when subjected to 
considerable -mag nification. Therefore, it is reasonable to assume that 
very near the leading edge (S— tO) the value of P would actually be 
greater than the calculated value of P at the given distance S. How-, 
ever, this shoiild have a negligible effect on the rest of the P curve ro 

and also on the total imp ing ement on the wedge, since the effect of a w 

stagnation point would be limited to a very small region about the leading 
edge. 


For S approaching very large values, the calculated values of P 
probably differ somewhat from act-ual values obtained in flight, because 
the analytical solution of the present report does not consider the 
existence of the bo-undary layer on the wedge surface. Since the botindary- 
layer thi.ckness increases -with the s-urface distance along the wedge, 
droplets which ing)inge at large distances from the leading edge act-ually 
■VTOTild have traveled throTJgh the boundary layer for some non- negligible 
interval of time. However, only a very small fraction of the total water 
droplets of a given size twinge -under this condition, and for large 
values of S the values of P are negligibly small. For example, in 
figure 4(a) such -would be the case for values of S greater than 8 or 
9 feet, where P < 0.0002 as conipared with Pq = 0.02. 

A preliminary stirvey disclosed a negligible effect of the free- 
stream static tempera-ture on the local impingement efficiency as a func- 
tion of the surface distance along the wedge (p against S). Values 
of P sud corresponding values of S were calculated for free-stream 
static teii 5 )era-tures of 4200, 440°, and 460® R, for a free-stream Mach 
number of 1.3, pressure altit-ude of 15,000 feet, droplet diameter of 
20 microns, and tangent of the semi-apex angle of 0.06. The values of 
P for the free-stream static tenperature of 420° and 460° R are within 
1 percent of the values of P at the free-stream static temperature of 
440° R. Since these calculations show that curves of P against S for 
the three values of free-stream static ten 5 >era-ture form practically a 
single c-urve when plotted to the -usual scales, no figures are presented 
herein to illustrate the effect of free-stream static ten^erature on 
iii 5 >lngement . Furthermore, the resiilts included herein, which eire calcu- 
lated for a free-stream static tenpera-bure of 440° R, may be used in the 
range of ten^era-ture from 420° to 460° R or possibly an even greater 
range of temperat-ure . 

The effect of the free-stream static pressure on the local impinge- 
ment efficiency as a f\mctlon of the surface distance along the wedge 
(P against S) is presented in figure 4(a) for pressijre altitudes of 
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sea level, 15,000, and 30,000 feet. Increasing the pressure altitude 
(decreasing the free-stream static pressure) increases the values of P 
at any distance S. For example, at S = 1.5 at sea level, p Is 0.0052j 
and at 30,000 feet, P is 0.0071. Since Pq (the value of P at 

S = 0) is equal to the sin a, where a is the semi-apex angle, the 
curves for the three pressvire altitudes have the same Tnn.xiTmTni value of p. 
Also, _ for the various pressure altitudes the extent of in^jingement along 
the wedge surface is essentially the same. 

The effect of the semi -apex angle o' of the wedge on the local 
ingjingement efficiency as a function of the surface distance along the 
wedge is presented in figure 4(h) for values of tan a from 0.02 to 0.10. 
Since the values of Pq are eqiml to sin a. Increasing the semi-apex 

angle of the wedge results in an increase of Pq. The surface extent of 
perceptible impingement (as characterized hy p « 0.01 Pq) does not 
vary as the wedge thickness is increased. 

The effect of free-stream Mach number on the p curve is presented 
in figure 4(c). The curves of P as a -function of S are given for 
free-stream Mach numbers of 1.2, 1.3, 1.4, 1.5, and 2.0. For the wedge 
semi-apex angle presented in the figure (tan a = 0.04), the value of 
M]_ = 1.2 is close to the shock- wave attachment Mach nimiber. The shock- 

wave attachment Mach number is a function of the wedge semi-apex angle 
and is defined as that Mach number below which the shock wave is detached 
from the wedge. An increase in the free-stream Mach number M- 2 _ results 

in an increased surface extent of perceptible in^pingement and, also in an 
incheaaed value of P at any given distance S (except at S = 0, where 
P = Pq = sin a and at S— >oo, where p — >0). The increase in the sur- 
face extent of perceptible in^ilngement is shown in figure 4(c), in which 
for free-stream Mach numbers of 1.2, 1.3, 1.4, 1.5, and -2.0 the surface 
extents of perceptible ingingement on the wedge (where p«0.01 Pq) are 
5.35, 6.05, 6.65, 7.20, and 9.4 feet, respectively. 

The effect of the droplet size on P is presented in figure 4(d) 
for droplet diameters of 10, 20, 30, 40, 50, and 100 microns. The surface 
extent of imrpingement and the values of ’P at any given distance S are 
considerably increased as the droplet size is increased. For exanple, 
for the semi-apex angle presented in the figure (tan 0 = 0.06), at 
S = 3 feet the values of P are 0.0000, 0.0033, 0.0106, 0.0182, 0.0244, 
and 0.0417 for values of droplet diameter of 10, 20, 30, 40, 50, and 
100 microns, respectively. The STirface extent of perceptible impingement 
has val ues of 1.5, 5.7, 11.2, 18.1, and. 2B.0 feet for droplet diameters 
of 10, 20, 30, 40, and 50 microns, respectively. As shown by the pre- 
ceding discussion and also by a comparison of the values of P as a 
function of S (fig. 4)., varying the droplet diameter from 20 to 
30 miicrons or from 30 to 40 microns is of the same order of magnitude in 
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its effect on the 3 against S curve as varying the pressure altitude 
from sea level to 30,000 feet or varying the free-stream Mach number from 
1.2 to 1.5. 

Total iiqpingement on wedge . - The effect of the free-stream Mach 
numher, the semi-apex angle of the wedge, the pressure altitude, and the 
droplet size on the total impingement on a wedge surface of infinite 
extent is given in figure 5. The total impingement is represented hy 
which is the ^ of the droplet having its trajectory tangent to the 

wedge surface (theoretically the tangent trajectory touches the wedge S 

sxirf ace at a point S — >oo ) . This can be obtained from the Integra- ^ 

tion 



or more directly from equation (l3) of the analytical solution, as 


^m ^ ^m ® 

where Xjj^ is given by equation (9) and is also defined in appendix A. 

The value for x^^ is obtained from the expression for x (eq. (?)) by 

allowing T— ><». Since the droplet which enters the flow field downstream 
of the shock wave at a distance ^ eqtial to (fig. l) theoretically 

in 5 )inges on the wedge surface only as the surface distance S approaches 
infinity, only droplets which have a value of ^ equal to or less than 
will impinge on the wedge surface of infinite extent. The rate of 

total water catch on one wedge sTirface in terms of is given as 

Wnx = 0.329^jj,ViW 

where is obtained in pounds per hour per foot of spein, is the 

free-stream velocity in miles per hour, and w is the liquid-water 
content in grams per cubic meter. Therefore, is directly proportional 

to the rate of total water catch on the entire wedge surface and is the 
rate of total water catch on one wedge surface per unit of span, free- 
stream velocity, and liquid- water content (in appropriate units). 

The variation of with the tangent of the wedge semi-apex angle, 

tan a, is shown in figure 5(a) for free-stream Mach numbers of 1.2, 1.3, 
1.4, and 2.0. As expected from the cTirves of 3 as a function of S, 
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the valtie of increases as tan a increases. However, the rate of 

the increase in with respect to tan 0 decreases with an increase 

in tan a. As can he seen from figure 5(a), increasing the free-stream 
Mach numher increases the value of for a constant value of tan a. 

The variation of with pressure altitude is presented in fig- 

ure 5(h) for two Mach numhers . In the range of pressure altitudes from 
sea level to 30,000 feet, the increase of with an increase in pres- 
sure altitude is approximately linear. The variation of with the 

droplet diameter d in microns is presented in figure 5(c). In the 
range of droplet diameter from 10 to 100 microns, as a function of d 

results in a curve that is very nearly a straight line when plotted on 
logarithmic paper. This linearity permits an accurate interpolation of 
when calculations have heen made for a few droplet diameters for a 

given value of wedge semi-apex angle, free-stream Mach numher, eind pres- 
sure altitude. 

Droplet velocities at impingement. - The variation of ,lm/^l 

(the ratio of the droplet inpingement velocity to the free-stream velocity) 
with the sTirface distance along the wedge is presented in figure 6 for 
three cases. The three cases given in figure 6 are representative of the 
resxilts when the drhplet diameter d is 20 niicrons and the pressure 
altitude is 15,000 or 30,000 feet. The curves of ixa^l ^ func- 
tion of S have characteristics s imi lar to the curves which present p 
€Ls a function of S. At S = 0, obviously, all the curves have in/^1 

equal to -unity; and as S is increased, the value of the velocity ratio 
rapidly decreases and asynq)totically approaches Vg/^l^ ratio of the 

air velocity downstream of the shock wave to the air velocity upstream of 
the shock wave. The top cuirve in figiure 6 ill-ustrates a typical situation 
for which very nearly unit y (tan 0 = 0.02 and = 2.0). 

The two lower curves (tan 0 = 0.06 and 0.10 at M^ = 1.3) are typical 

for cases where a stronger shock -wave produced hy a larger semi-apex angle 
results in decreased values of im/^l large values of S. 


Diamond Airfoil 

As sho-wn in the AMLISIS, the local impingement efficiency P at 
any point on the forward surfaces of a diamond airfoil (surfaces AB and 
AC shown in fig. 3) can he obtained directly from the results for the 
local impingement efficiency eis a function of the surface distance for 
wedges. The local inpingement efficiency p may he ohthined from fig- 
ure 4 or eq-uation (l7 ) . 
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In general, the results for the impingement on a diamond airfoil 
are presented in this report in terms of the total collection efficiency 
Ejj^ as a fimction of the scale parameter ilr, in the atten 5 >t to conform 

with the existing literature on the inpingement characteristics of air- 
foils. In the notation of the present report, the total collection effi- 
ciency Ejjj^ as stated in references 7 and 8 is defined as 


Em — 


^ui + Ull 


( 22 ) 


where T is the maximum thickness of the diamond airfoil, and | 

and are the absolute veQ.ues of the initial displacements from the 

leading edge (in a direction normal to free-stream direction) of the 
droplet trajectories which Impinge at the shotalder of the upper anii lower 
surfaces, respectively, of the diamond airfoil. For the symmetrical 
diamond airfoil at zero angle of attack, and will be equal. 

For the symmetrical diamond airfoil at an angle of attack, the tangent 
of which is equal to the thickness ratio, the value of is equal 

to zero. 


The scale parameter ijr is calculated for the diamond airfoil at 
supersonic speeds as for other airfoils at subsonic speeds. It is defined 
as 


♦ = (v) ( J) (2=) 

where c is the chord length of the diamond airfoil. The results pre- 
sented for inpingement on diamond airfoils lose essentially the same 
parameters as used for inpingement on wedges. The parameters are; free- 
stream Mach number, pressure altitude, diamond airfoil thickness ratio, 
and droplet size.' In addition, angle of attack and chord length are 
specified for the diamond airfoil. 

Total collection efficiency at zero angle of attack . - The variation 
of the total collection efficiency Ejjj with respect to the scale param- 
eter ijr is presented in figure 7 for the diamond airfoil at zero single 
of attack. The effects of pressiire altitude, airfoil thickness ratio, 
free-stream Mach number, and droplet size are shown. 

The effect of pressure altitude on the variation of Ejjj with ijr 

for a 2-percent-thick diamond airfoil is given in figure 7(a) for zero 
angle of attack, free-'stream Mach number of 1.4, droplet diameter of 
20 microns, and free-stream static temperature of 440° K. The pressure 
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altitudes presented in the figure are for sea level, 15,000 feet, and 

30.000 feet. The lines for constant values of chord length from 1 to 

' 20 feet are also included in the figure. For the diamond airfoil sub- 

jected to a constant-velocity supersonic flow field, the total collection 
efficiency Ej^ increases as the pressure altitude is increased when the 

chord length and the other variables are held constant. Considering Ejj^ 

as a function of scale parameter ijr where the droplet free-stream 
Reynolds number Re-i is held constant yields results similar to those 

DO 

for airfoils having rounded leading edges at subsonic speeds. With Re^ 
held constant, the veilue of decreases as \lr increases in the Tnnnnp-r 

Indicated in figure 7(a). As previoiasly mentioned, a con^jarison of the 
total collection efficiency and also the water impingement rate for a 
diamond airfoil at zero angle of attack with those for a MCA 0006-64 
airfoil at zero angle of attack is presented in appendix C. 

The effect of airfoil thickness ratio on the variation of E^ with 

is presented in figure 7(b). The curves are presented for 2- to 
14-percent-thick: diamond airfoils (in steps of 2 percent) at zero angle 
of attack, free-stream Mach nimiber of 1.4, pressure altitude of 

15.000 feet, droplet diameter of 20 microns, and free-stream static 
temperature of 440° R. The droplet free-stream Reynolds number is main- 
tained at a value of 453 for all the curves. The effect of increasing 
the airfoil thickness ratio is to decrease the total collection efficiency. 
The rate of decrease in the Ej^ with an increase in the airfoil thickness 

ratio becomes somewhat smaller as the airfoil thickness ratio increases. 

The effect of free-stream Mach number on the variation of total 
collection efficiency Em with scale parameter i]r is shown in fig- 
ure 7(c). The curves presented are for free-stream Mach numbers of 1.1, 

1.2, 1.3, 1.4, 1.5, and 2.0, 2 -per cent- thick diamond airfoil, zero angle 
of attack, pressure altitude of 15,000 feet, droplet diameter of 20 microns, 
and free-stream static ten^jerature of 440° R. The resTilts in the figure 
show that the collection efficiency increases as the frefe- stream Mach 
number Increases. ' However, the Increase in the total collection efficiency 
from the free-stream Mach number of 1.1 to 1.2 is considerably greater 
than the Increase in efficiency from a free-stream Mach number of 1.2 to 
1.3 eind from 1.3 to 1.4, and so forth. Calculations for 4- h.tk^ 6-percent- 
thlck diamond airfoils indicate the same trend. For the 2-percent-thick 
diamond airfoil the lowest Mach number presented in figure 7(c) (M = l.l) 
is quite close to the limiting Mach number for shock- wave attachment. 
Therefore, the rate of decrease of total collection efficiency with 
decreases in free-stream Mach number Increases as the Mach number 
approaches the shock-wave attachment Mach number as a limit. 
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The effects of droplet size and. chord length on the variation of 
total collection efficiency Ejj^ with scale parameter ijr are shown in 

figures 7(d) and (e) . Figure 7(d) presents curves of Em against ijr 
for a 2-percent- thick diamond airfoil, zero angle of attack, pressure 
altitride of 15,000 feet, free-stream Mach number of 1.4, and free-stream 
static tenperature of 440° R. The ciarves presented in the figure are for 
consteint values of droplet diameter (d = 10, 20, 30, 40, and 50 microns) 
eis well as for constant values of chord length (c =1, 2, 4, 8, and 
20 ft). Tncreetsing the droplet size greatly increases the total collec- 
tion efficiency. For exan^le, at c = 8 feet, the values of E^ are 

0.096, 0.310, 0.495, 0.625, and 0.711 at droplet diameters of 10, 20, 30, 
40, arid 50 microns, respectively. The rate of Increase in the total 
collection efficiency as the droplet diameter increases is less for the 
larger droplet sizes. This effect can also be observed in figure 7(e), 
in which curves of Em. as a function of ijr are presented for a 6-percent- 
thick diamond airfoil, zero angle of attack, pressure altitude of 
15,000 feet, free-stream Mach number of 1.3, free-stream static tempera- 
ture of 4400 R^ and for droplet diameters of 2, 10, 20, 30, 40, 50, and 
100 microns. 

A conparison of figures 7(b) and (c) (same droplet size and pressure 
altitude) shows that the effect on the total collection efficiency E^ 

of an increase in the free-stream Mach number from 1.1 to 2.0 is, in 
general, of the same order of magnitude as a decrease from a 14- to a 
2-percent- thick diamond airfoil. For example, for iJr of 1790 
(c = 8 ft) in figtire 7(b), the values of Em decreased from 0.310 to 
0.165 for an increase from a 2- to a 14-percent- thick diamond airfoil, 
respectively, a decrease of 0.145 in the value of E^^. For ilr of 1790 

(c c 8 ft) in figure 7(c), the values of Ej^ increased from 0.230 to 

0.400 for an increase in the free-stream Mach number from 1.1 to 2.0, 
respectively, an increase of 0.170 in the value of Ejjj. 

Another comparison of figures 7(c) and (d) (same airfoil thickness 
ratio and pressure altitiode) shows that the effect on the total collection 
efficiency E^ of an increase in the droplet diameter from 10 to 
50 microns is much greater than an increase in the free-stream Mach number 
from 1.1 to 2.0. For example, for c = 8 feet in figure 7(d) the value 
of Ejj^ increased from 0.095 to 0.710 for an increase in the droplet 

diameter from 10 to 50 microns, respectively, which is an increase of 
0.615 in the value of Em. As stated previoiosly, for c t= 8 feet 
(ilf = 1790) in figure 7(c) the increase in the value of Em is 0.170 for 
a corresponding increase in the value of the free-stream Mach number from 
1.1 to 2.0. For a constant value of the chord length, varyirg the pres- 
sure altitude has a re'latively small effect on the total collection effi- 
ciency (fig. 7(a)). 
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difl'mnnd airfoil . - The previous discussion is concerned with the totsil 
collection efficiencies for the diamond airfoil at zero angle of attack 
only. The effect of angle of attack on the total collection efficiency 
is presented in figure 8, in which the variation of Em with the angle 
of attack a is shown for 1- to 8-percent- thick diamond airfoils, inclu- 
sive, pressure altitude of 15,000 feet, free-stream Mach numher of 1.4, 
droplet diameter of 20 microns, and, chord lengths of 1, 2, 4, and 12 feet. 
The range in the angle of attack presented in f Igiire 8 is from a of 
zero to a of tan“^ (t/c). 

The toteil collection efficiency decreases slightly as the angle of 
attack increases. The slope of the curve of Ej^ against tan a is zero 

at a = 0, hy virtue of the symmetry of the diamond airfoil. As the 
angle of attack increases, the slope of the against tan a curve 

becomes negative, the rate of change becoming more pronounced as the 
angle of attack a increases. The effect of increasing the chord length 
of the diamond airfoil is to decrease slightly the rate of decrease of 
Ejj^ with respect to the angle of attack a. A significant point that is 

well- illustrated in figure 8 is that, for a constant value of chord length, 
there apparently exists an envelope of the family of Em against tan a 
curves that have the airfoil thickness ratio t/c as the parameter. This 
envelope ciurve presents, in terms of the angle of attack, the maximum 
total collection efficiency that can be obtained for a diamond airfoil of 
any thickness ratio, where the droplet size, free-stream Mach number, 
pressiire altitiide, and chord length are considered to be constant. 

For the diamond airfoil at a free-stream Mach number greater tha n 
the shock-wave attachment Mach number, the decrease in total collection 
efficiency with an increase in angle of attack is opposed to the trend 
experienced by roimded- leading- edge airfoils at subsonic Mach numbers 
(irrespective of the symmetry of the airfoil). For the latter type of 
airfoil at subsonic speeds the increase in total collection efficiency 
(as defined herein) with Increasing angle of attack is accounted for by 
the greatly increased impingement that occ-urs on the lower surfaces of 
these airfoils at angle of attack. The reduction of impingement occurring 
on the upper surfaces of these airfoils at angle of attack is more than 
balanced by the increased Impingement occurring on the lower surfaces. 

On the other hand, for a given diamond airfoil at supersonic speeds with 
an attached shock wave, the rate of increase with angle of atteick of the 
iD^iingement on the lower surface is less in magnitude than the rate of 
decrease with angle of attack of the impingement on the upper surface. 

This general trend for the diamond airfoil can be explained with the help 
of the figiire 5(a), in which the effect of tan a (tangent of the semi- 
apex angle) on the total impingement rate is presented for a wedge of 
infinite extent. 



22 


MCA TN 2971 


The forward upper surface of a diamond airfoil effectively acts as 
the finite portion of one surface of an infinite "weclge -which is decreasing 
its semi-apex angle (effectively the thiclhaess of -wedge) as the angle of 
attack of the diamond airfoil Increases. The forward lower surface of a 
diamond airfoil effectively acts as the finite portion of one surface of 
an infinite wedge that is increasing its semi-apex angle as the angle of 
attack of the diamond airfoil increases. In figure 5(a) the increase in 
which is exactly proportional to the total impingement rate (see 

appendix A), hecomes smaller -with an increase in tan a (i.e., 
d2 ^jjj/d(tan a)2 < o) for all possible values of the semi -apex eingle cr. 
For example, in figure 5(a) for the curve of = 1.4, the impingement 

on an 8-percent-thick diamond airfoil of large chord extent and at zero 
angle of attack -would he very nearly proportional to -twice the value of 
at -tan ci of 0.08 (2^j^ = 0.146 ft). The impingement for the same 

airfoil under -the sarnie cond i tions and at an angle of attack of 2.291° 

(tan a = 0.b4) -would he very nearly proportional, to the sum of at 

tan a of 0.04 and 0.12 (^m,u + ^m,2 = 0-046 + 0.0905 = 0.1365 ft). 

Total collection efficiency as function of modified inertia param- 
eter . - The total collection efficiency is presented in figure 9 as 

a function of the imcdified Inertia parameter Kq, as suggested in refer- 
ence 14. The relation for Kq is 

Ko = (XAs)K 

where K is -the Inertia parameter and is defined as 

P 

and A is the maximum distance of travel of a droplet projected into 
still air with the free-stream velocity Vp. The term Ag is the value 

of the maximum distance of -travel A -rfien Stokes' law is assumed for the 
drag force on the droplet. Since K is equal to Ag/c (ref. 4), Kq 

may also he -written as 

Kq = (a/As)K = A/c 

The free-stream droplet Reynolds numibers used herein in calculating 
Kq range -up to 2104. The empirical drag law (eq. (5)) \ised in this 

report for the droplet motion in the mso-vliig reference frame is valid up 
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to the Reynolds number of approximately 500. Therefore, for the values 
of Kq presented in figure 9 (and also fig. 10), the values of (A/\g) 

■were obtained from table I of reference 4. 

The correlation of the total collection efficiency as a function of 
the modified inertia parameter Kq is shown in figure 9(a) for pressure 

altlt'udes of sea level, 15,000, and 30,000 feet, and for droplet diameters 
varying from 2 to 100 microns. For the 6-percent-thick diamon,(l airfoil 
at M-j_ = 1.3 and zero angle of attack presented in figure 9(a), the 

droplet free- stream Reynolds number Re^ varies from 42.1 to 2104. The 
values of Ej^ for a given thickness ratio, Mach number, and angle of 

attack generally form the baais for a single curve with a small amount of 
scatter existing in the higher range of value of Kq. In the lower range 

of Kq, for all the values of Re^^ the plotted points have negligible 

scatter. The small scatter observed is possibly due to the existence of 
a veiy slight trend of the Em against Kq curves -with Rei for the 
higher values of Kq. 

The effect of increasing the thickness ratio is presented in fig- 
lore 9(b) for 2-, 6-, and 12-percent thick diamond airfoils. Increasing 
the thickness ratio displaces the Ejj^ against Kq curve toward larger 

values of Kq. Changing the thickness ratio of the diamond airfoil does 
not alter the shape of the c'urve itself. 

The effect of the free-stream Mach number on the variation of Em 
■with Kq is presented in figure 9(c) for the 2-percent- thick i nmnuf^ 

airfoil, droplet diameter of 20 microns, pressure altitude of 15,000 feet, 
and for free-stream Mach nxmibers of 1.1, 1.2, 1.5, anfl. 2.0. Increasing 
the Mach number displaces the entire curve of Em against Kq to'wards 
smaller values of Kq. As the value of M^ increases, the rate of dis- 
placement of the c^urve -with Increasing Mp becomes smaller. The dis- 
placement of the curve obtained by increasing the Mach number from 1.1 
to 1.2 is more than that obtained by increasing the free-stream Mach 
number from 1.2 to 1.5 and from 1.5 to 2.0. 

The total collection efficiency for a diamond airfoil can also be 
considered as a fimction of another nidified inertia parameter defined as 
a relative modified inertia parameter. The relative modified inertia 
parameter Fq is based on the velocity change of the air across the shock 

wave rather than on the free-stream air velocity. This modified inertia 
parameter and its effect on the total collection efficiency of a diamond 
airfoil are presented and disciissed in appendix D. 
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SUMMART OF RESULTS 

This report presents an analysis of the prohlem of impingement of 
water droplets on a wedge and a diamond, airfoil at supersonic speeds and 
with attached shock waves. When a suitable eng)irical relation is used 
for the drag coefficient of a sphere, the analysis allows a closed-form 
integration of the eqmtions of motion for the water droplets. The 
integration results in analytical expressions for the equation of the 
trajectories, the droplet velocity at aiay point on the trajectories, the 
local impingement efficiencies, and the total rate of impingement. The 
results of the calculations of the rate, extent, and distribution of the 
impingement of water droplets on wedges and diamond airfoils are summa- 
rized briefly as follows (the Mach number referred to is the free-stream 
Mach nimmber, which is greater than the attachment Mach number for the 
wedge or the diamond airfoil): 

1. At a given value of Mach number, droplet size, and pressure alti- 
t\ade, the local impingement efficiency as a function of the dimensional 
surface distance is the same for both the wedge and the d ^ a-mnnd airfoil 
at zero angle of attack, provided the tangent of the semi-apex angle of 
the wedge is equal to the diamond airfoil thickness ratio. 

2. For any Mach number, pressTore altitude, and droplet diameter, the 
value of Po is equal to the sine of the semi- apex angle for wedge or 
diamond airfoil. (Pq is the value (maximum) of local impingemient 
efficiency as distance from leading edge to point of impingement 
approaches zero.) 

3. The effect of the free-stream static temperature on the local 
impingement efficiency and total collection efficiency is negligible for 
temperatures from 42CP to 460° R. 

4. At constant values of Ifech numiber, droplet size, and semi-apex 
angle of the wedge or corresponding thickness ratio of the diamond air- 
foil, an increase in the pressure altitude increases slightly the local 
impingement efficiencies and total collection rates on wedges and diamond 
airfoils; but an increase in the pressure altitude has a negligible 
effect on the surface extent of perceptible impingement. 

5. At constant values of lyfeich number, droplet size, and pressure 
altitTode, increasing the thickness ratio of the diamond airfoil or corre- 
sponding semi-apex angle of the wedge increases the local impingement 
efficiency, has a negligible effect on the stirface extent of perceptible 
Impingement, and decreases the total collection efficiency of the diamond 
airfoil. 

6. At constant values of droplet size, pressure altitude, and semi- 
apex angle of the wedge or thickness ratio of the diamond airfoil, an 
increase in Mach numober increases both the STorface extent of impingement 
and the value of the local Imroingement ' efficiency. 
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7. At constant values of pressure altitucLe, semi-apex angle of the 
wedge or thickness ratio of the diamond airfoil, and Mach number, an 
increase in the droplet size increases considerably the stirface extent 
of perceptible impingement, the local inpingement efficiency, and total 
impingement rate. 

8. The variation of total collection efficiency of the diamond air- 
foil at zero angle of attack as a function of the scale peirameter for 
constant values of the droplet free-stream Reynolds number is similar in 
form to that for subsonic airfoils. 

9. The total collection efficiency of the diamond airfoil decreases 
slightly as the angle of attack increases. 

10. For a diamond airfoil of a given thickness ratio and Mach 
number, the values of total collection efficiency for a wide raiage of 
values of droplet free-stream Reynolds number coinprlse a single curve 
when plotted against the modified inertia parameter. The effect of 
increasing the thickness ratio or decreasing the Mach number is to dis- 
place the entire collection efficiency against modified inertia parameter 
cTnrve in the direction of larger values of modified Inertia parameter. 


Lewis Flight Propiilsion Laboratory 

National Advisory Conmittee for Aeronautics 
Cleveland, Ohio, April 21, 1953 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 
droplet radius, ft 
airfoil chord length, ft 
drag coefficient, d im ensionless 

drag force on spherical water droplet, Ih 
droplet diameter, microns 

total collection efficiency (defined hy eq. (22)), dimension- 
less 

p a% X 

relative inertia parameter, (s/9) ^ ^ ^ =, — dimensionless 
relative modified inertia parameter, ~ ^ 

acceleration due to gravity, 32.2 ft/sec2 

/ / X Pw*^^l ^s 

inertia parameter, (2/9) ■ — = — , dimensionless 

^ C 


KO 

Ml 

m 

o^l^ 112,113 

5l 

R 

Re^ 


modified inertia parameter, (^/^g)K =; — , dimensionless 


free-stream Mach number 

en5)lrical constant (used in eq. (5) ), 2/S 

constants of flow field (defined hy eqs. (I8a), (I8h), and (I80), 
respectively) 

free-stream static pressure, Ih/sq ft 
gan constant, 53.3 ft-lh/(lb)(c>F) 

droplet Reynolds numher relative to air behind shock wave, 
2ap2U/p2, dimensionless 
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Re 2 _ free-stream droplet Reynolds number, 2a.p^2./\^±) dimensionless 

S distance to point of impingement measured from leading edge 

for a watCT droplet whiela enters the flow field behind the 
shock wave at a distance ^ above the leading edge 
(eg. (l5) and fig. l), ft 

T maximimi thickness of airfoil, ft 

T/c ratio of maximum thickness to chord length of airfoil, dimen- 

sionless 

t time, sec 

tj free-stream static temperature, OR 

U magnitude of droplet velocity relative to air velocity down- 

stream of shock wave, jv^ - V 2 I, fps 

droplet velocity, fps 

VjL free-stream velocity, fps or mph 

V 2 air velocity downstream of shock wave, fps 

W rate of water catch, lb/(hr)(ft span) 

Wp local rate of water catch, lb/(hr)(sq ft surface) 

w liquid- water content of air, g/cu m ' 


X 

^,s 

a 



U dt, displacement of droplet in moving frame of reference 


(relative to air stream), where t = 0 the instant the water 
droplet intercepts the shock wave 


maximum value of displacement x as t— as a limit 

(maximum distance of travel of droplet projected into still 
air with relative velocity U^), ft 


value of maximum distance of travel obtained by assuming 

. Stokes’ law for drag force on droplet, ft 

angle of attack of airfoil, deg 
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local impingement efficiency, d^/dS, dimensionless 

maximum value of P as S — >0; Pq = sin a, dimensionless 

ratio of specific heats (1.4 for air) 

empirical constant (used in eq. (5)), 0.158, dimensionless 

initial displacement of droplet trajectories from leading 
edge in direction normal to free-stream direction (eq. (13)) 

maximum value of initial displacement ^ for a droplet 
obtained when its trajectory is tangent to wedge surface 
(theoretically &e S— ><») 

absolute values of initial displacement from leading edge 
{in direction normal to free-stream direction) of droplet 
trajectories which impinge at shoulder of upper and lower 
STorfaces, respectively, of diamond airfoil 

distance along shock wave measured from wedge apex to point 
where droplet trajectory Intercepts shock wave 

shock-wave augle 

angle between free-stream velocity vector and droplet 

velocity vector 

maximum distance of travel of droplet projected into still air 
with free-stream velocity Vj, ft 

value of distance of travel A obtained by assuming Stokes' 
law for drag force on water droplet, ft 

dynamic viscosity at free-stream static ten 5 )erature, 
(lb)(sec)/sq ft 

dynamic viscosity at static tengjerature downstream of shock 
wave, (lb)(sec)/sq ft 

angle between free-stream velocity vector and initial 

relative velocity vector (defined by eq. (13a)) 

components of droplet displacement referred to wedge surface 
(defined by eq. (12)), ft 
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water-droplet mass density, 1.9398 slugs/cu ft 
free-stream mass density of air, slug/cu ft 
Pg mass density of air downstream of shock wave, slug/cu ft 

a semi-apex 6ingle 

a' angle between droplet velocity vector and velocity 

vector V2 

T dimensionless time variable, 

cp phase angle, tan“^Rej.^^“V2g“^/2)^ 0< cp< jt/2 (eq. (9)) 

i|f scale parameter, (9c/a)(pj^/p^) 

ratio of air velocity downstream of shock wave to free-stream 
velocity 

<0 ratio of initial droplet relative velocity to free-stream 

velocity U^/V^ 

Subscripts: 

i initial (at shock wave) 

im impingement 

Z lower 

m. maximum 

u upper 

Barred symbols denote vectorial qxiantities. 
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APPENDIX B 


CLOSED-EORM INTEGRATION OF EQUATION OF MOTION 
(RELATIVE TO AIR FLOW DOWNSTREAM OF SHOCK WAVE) 

This appendix provides the varioTis steps necessary to the closed- 
form integration of equation (6), which is the differential equation of 
motion of the water droplet in the moving reference frame. The differ- 
ential equation of motion (6) can he rewritten as 


dt 



(1 + cAiVa“) 


where for convenience 


(Bl) 


ROj, = 


Aj^Ua 


BP2 

1^2 


3 ^2 


^ = 8 ~ 
w 


(B2) 

(B3) 


(B4) 


(a complete list of symbols is given in appendix A. ) Algebraic simplifi- 
cation and rewriting result in the form 


A., dU 

dt ; 

U(1 + A^U®) 

where 


(B5) 


Ag = (Ai/A2),(a2/24) ' (B6) 


and 


A4 = eAj^V 


(B7) 
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Equation (B5) 1b not readily Integrable in its present form. By letting 

p = U® (B8) 

and 


equation (B5) becomes 


P’ = p/Pi = 


(B9) 


dt = - (A3/m)(dp*/p’) + (-As/m) 


A4Pi dp' 

(1 + 


(BIO) 


This form of the equation is readily integrated, and upon resubstituting 
the relations (B8) and (B9), there is obtained the following expression 
for the velocity of the droplet as a function of the time t: 


U = 



(Bll) 


where B^^ is the constant of the Integration. 

Since U = equation (Bll) can be rewritten in integral form as 


E 

where Bg is the second Integration constant. Equation (B12) requires 

siDipllflcation before a closed-form integration can be performed. Con- 
sider the following substitutions: 


- “l%] “ + ®2 


(B12) 




(B13) 

A_ „ -l/m A^ -l/m 

*5 = ^ “i (VAi) = V 

(B14) 

^ 1 

Ag = (U^ A^)” expC-Bj^m/Aj) 

(B15) 
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By visirig the Substitutions given By (E13), (B14), and (B15), it is 
possible to "write eqmtion (B12) as 

Jdx = {AqB^ - <3y + Bg (B16) 

Before further steps can he taken in the closed-fomi integration of the 
equation of motion of the droplet, it is necessary to detennine the value 
of the eii 5 )lrical constant m. It is noted that in approximatlDg the 
curve of the drag coefficient as a f"unction of the local Reynolds number 
by a relation of the form given by 


Cj) = (24/Rej,)(l + eRe^“) (S) 

it is possible to consider that "tdie value of m is 2/3 and the value of 
e is 0.158. That the approximation of the experimental curve by the 
enplrical relation is very good in the range of Reynolds numbers from 
about 0.5 to 500 can be seen by referring to figure 2, which presents a 
graph of the eng)irical relation along "with the drag- coefficient data of 
references 4 and 12. 

The "Use of m = 2/3 in equation (B16) along "with the substitution 

e^ = q (B17) 

and the tise of formulas of integration given on pages 16 and 17 of refer- 
ence 15 allow eq"uatlon (B16) to be integrated as shown in the following 
steps; 



Ag dq 
q(Agq - 1) 


3/2 ®2 



A 3 dq 


q(Agq - 1 ) 




A 3 dq 


(Agq - 1) 




H 

- Pfl - = - - tan-lVAgl - 1 (B18) 

2% - 1 
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If the original time or independent variable t is reintroduced into 
equation (B18) and the substitutions are made for Ag and Ag, the equa- 
tion has the form 


(A^)3/2 (Bg-x) 
3A^ 


(^^^^2/3)l/2 

I f2(t-\)] 


tan' 


-1 





(B19) 


Substituting for Ag and A^ (except in the exponent of e) in equa- 
tion (B19) results in 



Equation (B20) is the integrated equation with undetermined integration 
constants for the motion of the water droplets relative to the air veloc- 
ity behind the shock wave. The integration constants are determined from 
the botondary conditions^ which are 


U = Ui 


at t = 0 


X = 0 

The substitution of the boundary conditions, and thus the determination 
of the integration constants B^_ and Bg, resTolts in the final form of 

the integrated equation of motion for the water droplets' relative to the 
air velocity downstream of the shock wave as follows (in the frame of 
reference moving with the constant velocity Vg): 
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viler e 

T dimensionless time variable, (3ii2/pv a^)t 

initial value of local relative Reynolds number, Re^, EapgUj^/p^ 

a droplet radius, ft 

vater density, 1.9398 slugs/cu ft 

pg density of air behind shock wave, slug/cu ft 

c en^jirlcal constant of relation of drag coefficient as function 

of Reynolds number, 0.158 

The final form of the corresponding equation for the relative velocity 
of the water droplet as a function of the dimensionless time variable 1 
is obtained from equation (BU) as 



It can be seen from equations (?) and (8) that as t approaches infinity 
the value of U approaches zero and that a limit exists for the value of 
X as t approaches infinity. This limiting value of x is 

where 

cp = tan-1 0 ^ | 
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APPENDIX C 


COMPARISON OF COLLECTION EFFICIENCY AND IMPINGEMENT BATE FOR 
DIAMOND AIRFOIL raTH THOSE FOR NACA 0006-64 AIRFOIL 

A comparison of the total collection efficiency as a function of a 
modified inertia parameter Kq for a diamond airfoil at supersonic 

speeds (attached shock wave) with that for the NACA 0006-64 airfoil 
(ref. 9) at free-stream Mach numbers less than critical is presented in 
figmre 10. Both airfoils are 6-percent thick and are at zero single of 
attack. It is in^jortant to keep in mind that figure 10 does not allow 
a coirparison of the two airfoils at the same Mach number, eind such a 
comparison cannot be made, because the analysis for the NACA 0006-64 
airfoil is not valid above the critical Mach number and that for the 
diamond airfoil is not valid below the shock- wave attachment Mach number. 

It is, however, of value to consider a conparison of the two airfoils 
with each operating within its appropriate speed range. In figure 10 
the rate of increase in Ejj^ with respect to Kq is greater for the 

diamond airfoil than for the NACA 0006-64. For a constant droplet size, 
presstire altitude and tenperatinre, and free-stream Mach number resulting 
in a constant value of A, Kq varies inversely as the chord length 

(Kq = A/c). In general, therefore, the rate of decrease in E^ with 

increasing chord length is greater in magnitude for the diamond airfoil 
than for the NACA 0006-64 airfoil. 

A conparison of the rate of total water catch per u ni t span for the 
diamond airfoil with that for the NACA 0006-64 airfoil can be obtained 
if values (necessarily different) for the free-stream Mach numbers for 
the two airfoils are assigned. AssTxme that both airfoils are of the same 
thickness ratio and chord length. Further, assimae thAt the diamond air- 
foil and NACA 0006-64 airfoil encounter identical icing conditions of 
droplet diameter (20 microns), pressure altitude and temperature 
(15,000 ft and 440° R), and liquid-water content (0.5 g/cu m), but have 
free-stream Mach numbers of 1.5 and 0.75, respectively. For the static 
tenperature eissumed, the Mach numbers of 1.5 and 0.75 correspond to speeds 
of 1051 miles per hour and 526 miles per hour, respectively. Therefore, 
the magnitudes of the inertia parameter K and the free-stream droplet 
Reynolds ntmiber Re^^, for the diamond airfoil are twice as great as for the 
NACA 0006-64 airfoil. Varying the chord length from 1 to 20 feet produces 
a change in the value of Kq from 0.-386 to 0.0193 for the diamond air- 
foil and from 0.266 to 0.0130 for the NACA 0006-64 airfoil. This vari- 
ation in Kq for both airfoils results in values of Em which are of 
the same order of magnitude. For the given icing conditlqns, the follow- 
ing table lists for chord lengths of 4 and 20 feet the various pertinent 
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parameters and variables ^ including the rate of total water catch on the 
airfoil per unit span for both the 6-percent-thick diamond aifroil and 
the HACA 0006-64 airfoil; 


6 

-Percent-thick diamond 

airfoil at M^ = I 

L.5 (Vi = 1051 nph) 

Chord 

Free-stream 

Inertia 

X 

Ifodified 

Total 

Rate of water 

length. 

droplet 

parameter. 


inertia 

collection 

catch, W, 

c, ft 

Reynolds 

K 


parameter. 

efficiency. 

lb/(hr) 


number. 



Ko 

Em 

(unit span) 


Rei* 






4 

485 

0.542 

0.178 

0.0965 

0.432 

17.9 

20 

485 

.108 

.178 

.0193 

.111 

23.0 


KACA 0006-64 at M^ = 

0.75 

= 526 nph) 


4 

243 

0.271 

0.246 

0.0665 

0.330 

6.85 

20 

243 

.0542 

.246 

.0130 

.128 

13.3 


* Values of viscosity of air obtained from linear variation with tenpera- 
ture given in reference 16. 


The rate of total v&ter catch on the airfoil per unit span is calculated 
from 


Wjj^ = 0.329 EjjjTV2W 

where T is the TnaximuTn thickness of the airfoil in feet, w is the 
liquid- water content in grams per cubic meter, and Vj_ is the free- stream 

velocity in miles, per hour. For this particular example, the table shows 
that for the 4-foot chord, the d i amnruT. airfoil has a somewhat larger 
value of Ejj^ than the KACA 0006-64 airfoil; and for the 20-foot chord, 

the opposite is the case. Conparison of the two airfoils for a given 
chord length shows that, as expected, the effect of the total collection 
efficiency on the rate of water catch is small when conpared irtth the 
free-stream velocity ratio chosen (2;l). The most significant comparison 
to be obtained from the table is +-hat, for the diamond airfoil, increasing 
the chord length by a factor of 5 (from 4 to 20 ft) results in an increase 
of only 28 percent in the rate of water catch; whereas, for the WACA air- 
foil, a like change in the chord length resiilts in an increase of 94 per- 
cent. This difference is the result of the fact that, unlike the local 
impingement on the rounded leading-edge airfoils at free-stream Mach 
numbers less than critical, the local impingement at a given point on the 
surface of fliflTnnnfl airfoils does not vary with the chord length of the 
diamond airfoil, which is at a supersonic Mach nirnnher above the shock- 
wave attachment Mach number. 
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APPENDIX D 


TOTAL COLLECTION EFFICIENCy OF DIAMOND AIRFOILS AS FUNCTION 
OF RELATIVE MODIFIED INERTIA PARAMETER Fq 

For air-flow fields that contain a shock wave, such as those con- 
sidered herein, the water droplets upon crossing the shock wave stiddenly 
have a velocity relative to the air. For the diamond airfoil with an 
attached shock wave, this initial relative velocity is the same for all 
the droplets entering the air-flow field downstream of the shock wave. 

As discussed in previous sections, this common initial velocity may 

he considered to he the initial velocity of projection of a droplet in 
a reference frame having no air motion in it and moving relative to the 
fixed frame of reference at the constant velocity of the air downstream 
of the shock wave. Therefore, it is of interest to define a new set of 
inertia parameters, F and Fq, based on the motion of the droplet in 

this moving reference frame. In the expressions for Re^^, K, and Kq, 
if the initial relative droplet air velocity is substituted for Vj_, 

then 


p U 2a 
Rer,i = — 


F = K (Ui/Vi) 


(2/9) 


P a^U. 
*^w i 


X 

m,s 

c 


and 



where Xjjj (eq. (9)) is the Tnfl.x-iTtnm distance of travel of a droplet idiich 
is projected into still air with the relative velocity Ui, and Xm^s 
is that value of the maximum distemce of travel when Stokes ' law is 
assinned for the drag force on the droplet. 


The total collection efficiency Ejjj is presented in figinre 11 as a 
function of the relative modified inertia parameter Fq. The curves of 
Em against Fq have exactly the same form as the curves of Em against 
Kq, except that the use of Fq resiilts in a displacement of the entire 
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curve towards STna.13.er numerical values of the abscissa. The effect of 
diamond airfoil thiclmess ratio and of free- stream Mach numher on Em 
as a function of the relative modified inertia parameter Fq is presented 
in figures ll(a) and (h), respectively. The displacement of the Ejj^ 
against Fq curves obtained by varying the thiclmess ratio (fig. 11(a)) 

or by varying the free-stream Mach number (fig. 11(b)) is greater than 
the displacement of the Ej^ against Kq curves obtained by the same 

procedures (figs. 9(b) and (c)). 
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M 



Figure 2. - Comparison of eng)irical relation for drag coefficient as 
function of droplet Reynolds nuiiber vlth experimental values. 
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= 1 + 0.158(Hej,)^'^^. 
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1 ^ 

CO 


Shook Eipaiielon B^panBlon 



^a) Ai35le of a+rtack; a < ‘tan "^(T/o) , (h) Angle of attack, ct - tan"^(T/o). (o) Angle of attack, a > tan"^(T/o) . 


Figure 5. - Gohemtlc (ilagran of dlanond. airfoil in BuperBocio flow vith attached shook varea and at angle of attack. 
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DlBtanoe along surfac# of wedge (frooi leading edge), 3, ft 

(a) Sffeot of presBure altitude. Droplet diameter, 20 mlorons; freo-etream Maoh nuater, l.ij tangent of eeml-apax 
angle, 0.02. 

Figure i. - Local implngeoent efflolena^ on wedge as funotlon of dlatanoe along surface of wedge. Pree-stream static tee 

ruKT«B^in4k . AjtnO P 


Local Impingement efficiency 
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0 2 4 6 8 10 12 
• Distance along siirface of wedge (from leading edge), S, ft 


(b) Effect of wedge semi-apex angle. Droplet diameter, 

20 micronsj free-stream Mach number, 1.3; pressure 
ailtltude, 15,000 feet. 

Figure - Continued. Local Impingement efficiency on wedge 
as function of distance along siu?face of wedge. Free-stream 
static temperature, 440® R. 
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Tangent of aeml-apejc angle, tan o 

(a) Effect of tangent of seml-^pex angle. PreaHure altitude, 15,000 feet> droplet diameter, 

20 mlcrona. 

Figure 5. - Total ia^lngement rate V)^ for vedge of infinite extant, Pree-atream atatic 'ten 5 >era 
ture, 440° E. o 0.329 $_VtX. 
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Total coUactloa efficiency. 


oi 

CO 



(1)) affect of thlckaesB ratio of ilaaond airfoil*. Droplet dlamerter, 20 mlcronfl) frea-stream Mach mpiher, 
1.4) free-fltreaiB droplet Raynolda mniber, 463) ptreBsure altitude, 15,000 feet. 

nguro 7. - Continued, Total oolleotlon efficiency of diaBond alrfollo ae function of Boale parameter. Ifoe- 
atream atatic temperature, 440® B) angle of attaok, 0®. 


1 


29A2 


MCA !ffl 2971 





BoaJ.e pfiraaeiar, ^ 


(o) Effacrt of frea-straa* Modi number. Diamond airfoil thlckneoBi 2 percent) droplet 
dlameterj 20 microns) pressure altitude, 16,000 feet. 

figure 7. - Continued, total collection efficiency of diamond nirfoilfl aa function of seals 
porematar. rree-atroam static tetuperatura, 410° R) angle of attack, 0°, 





8 1 2 -t 6 8 10 20 40 60 80 lOOxlO^ 

floalB paraneter, 

(d) Bffect or droBlet site chord length. Dlajwsnd alxfoU thlalmaaa, 2 percent; frae-stcreajn Mach mmiber, 1.4; preBsure alti- 
tude, 15,000 feet. 

Figure 7. - Continued. Total collection efficiency of dlaaond alrfollfl an function of scale parameter. Froo-etream atatlo teiw 
perattrre, 440® Hj angle of attack, 0^. 
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Figure 8. - Variation of total collection efficiency vlthi angle of attack. Free-otream . 
static teB 5 >erat\jre, 440° Bj droplet dlanffiter, 20 microns} free-stream Mach number, 1.4} 
pressxrre altitude, 15,000 feet. 
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Droplet dieaeter. 
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Modified Inartla paraiaeter, Kq 
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(a) Correlation of preesnre altitude and droplet alse. Diamonii airfoil thlctaieaa, 6 percent; freo-etream Ifeoh 
nuoiiarj 1.3. 

Plgore 9. - Total Golloctlon eTficlaucy of dlaaond airfoils as fanctlon of Modified Inartla paranoter* Frea-stream static 
teim>eratiire| 440° ang.1.fl of attack^ 0^* 






Total collection efficiency, 
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Modified Inertia parameter, Kq 


(■fa) Effect of diamond airfoil thickness ratio. Eree-stream Mach number, 1.4j 
droplet diameter, 20 mloronsj pressure altitude, 15,000 feet. 

Figure 9. - Continued. Total collection efficiency of diamond airfoils as func- 
tion of modified inertia parameter. Free-stream static temperature, 440° E; 
angle of attack, 0°. 
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Total coiLlection efficiency 



.001 .002 .004 .006 .01 .02 ,04 .06 .1 .2 .4 ,6 1 2 

Modified inertia parameter^ Kq 


Figure 10. - Con^parison of total collection efficiency aa function of modified 
inertia parameter for diamond airfoil at superaonic speeds (attached shock ^, 76 ) 
■with that for HACA 0006-64 airfoil at frea-atream Mach numbers less than criti- 
cal. Airfoil thiclmeBs, 6 percentf angle of attack, 0°. 
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Belatlve molified inertia parameter, Pq 

(a) Effect of diamond airfoil thickness ratio. Free-etareeun Mach number, 1.4, 

Figure 11. ' Total collection efficiency of diamond airfoils as function of relative modified inertia 
parameter, Free-atream static teii^perature, 440° E> angle of attack, 0°. 
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Figure 11. - Oonoluded. Total ooUeotlon effloienoy of diamond airfoils as fonotion of 
relative modified inertia parameter. Free-atream static temperature, 440° R; angle 
of attack, 0°. 





